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We present a consistent calculation of the structure functions within a light-front 
constituent quark model of the nucleon. Relativistic effects and the relevance of 
the covariance constraints are analyzed for polarized parton distributions. Various 
models, which differ in their gluonic structure at the hadronic scale, are investi- 
gated. The results of the full covariant calculation are compared with those of 
a non-relativistic approximation to show the structure and magnitude of the dif- 
ferences. It is also shown how measurements of transversity in doubly polarized 
Drell-Yan lepton pair production are a clearcut sign of covariance requirements for 
the spin. 



1 Introduction 

We report here on calculations of the polarized structure functions within a 
covariant quark model of the nucleon. The framework used is based on a quark- 
parton picture we have developed to describe the structure on the nucleons in 
the deep inelastic regime El. The central assumption in this procedure □ is 
the existence of a scale, ^Lq, where the short range (perturbative) part of the 
interaction is negligible, therefore the glue and sea are suppressed, and the 
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long range (confining) part of the interaction produces a proton composed of 
three (valence) quarks only. Jaffe and RossH proposed thereafter, to ascribe 
the quark model calculations of matrix elements to that hadronic scale /Iq. 
For larger Q 2 their Wilson coefficients will give the evolution as dictated by 
perturbative QCD. In this way quark models, summarizing a great deal of 
hadronic properties, may substitute ad hoc low-energy parameterizations. 

However, it has been quite evident, since the original formulation of many 
of these models, that relativistic effects to the nucleon wave function, as well as 
covariance requirements, are needed even for a phenomcnological description of 
the structure of hadrons. To accomplish with this aim we develop a constituent 
quark model in the light-front realization of the Hamiltonian dynamics whcxp 
covariance can be incorporated in a rather transparent and elegant manner 

We apply this formalism to the calculation of polarized parton distribution. 
The role of initial soft gluons in the description of the parton distribution (and 
especially gluon polarization) is analyzed. Finally we show how differences 
between transverse and longitudinal polarization are a direct consequence of a 
light-front treatment of the spin and they may be detected experimentally in 
a rather model independent way. 

2 Polarized parton distributions at the hadronic scale 

The parton distributions at the hadronic scale /1q are assumed to be valence 
quarks and gluons. We will introduce the gluons later on in a phenomcnological 
way, let us therefore now restrict our discussion to the qiiarks. The quark 
distribution is determined by the quark momentum density EJ, 

where k + = ^/k 2 + to 2 + k z is the light-cone momentum of the struck parton, 

M/v and m q are the nucleon and (constituent) quark masses respectively and 
nj(k 2 ), rt^(k 2 ) represent the density momentum distributions of the valence 
quark of g-flavor aligned (|) or anti-aligned Q) to the total spin of the parent 
nucleon. 

nl$(k 2 ) = (N,J Z = +l/2\J2 1 + ( 2 ' )Tf 1 + ( 2 ' K % (k^k t )|iV,J 2 =+l/2). 

i=l 

The distributions (0) have been evaluated in the past making use of non- 
relativistic constituent quark models while, in the present investigation, we 
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want to improve their description by including relativistic effects as dictated 
by a light-front formulation of a three-body interacting system. n 

We just outline here the basic features of this formulation (seecl for a re- 
view). In light-front dynamics the intrinsic momenta of the constituent quarks 
(fej) can be obtained from the corresponding momenta in a generic ref- 
erence frame through a light-front boost (ki = CJ 1 (Ptot) Pi, Ptot = Si=i Pi) 
such that the Wigner rotations reduce to identities. The spin and spatial 
degrees of freedom are described by the wave function: 

* = -^«(P-p)x(ki,A*i,...,k8,A*3) (3) 

where /i, refer to the eigenvalue of the light-front spin, so that the spin part of 
the wave function is transformed by the tensor product of three independent 
Melosh rotations: i?^(kj, rm) U, namely — n® i=i ^uO s i' m i)- Finally) 
the internal wave function is an eigenstate of the baryon mass operator M = 
M Q +V, with M — Yli=i — Eti \/kf + to?, where the interaction term V 
must be independent on the total momentum P and invariant under rotations. 

In the present work we will discuss results of a confining mass equation of 
the following kind 



(M + V) Vo,o(£) = \ 1^ \/k? + to? - j + «,e I Vo,o(0 = MVo,o(0 , (4) 




where kj = 0, £ = y p 2 + X 2 is the radius of the hyper-sphere in six 
dimension and p and A are the intrinsic Jacobi coordinates p = (ri — r^j/V^i 
A = (ri +r2 — 2 ra)/ \/6- The choice of the mass operator (^) has been motivated 
by the fact that it combines a simple form with a reasonable description of the 
baryonic spectrum LrEL 

The intrinsic wave function (disregarding the color part) of the nucleon 
can be written: 



\N, J, J„ - +1/2) = Vo,o (0 y fo,o,o] [XAfs^MS + Xma<Pma] I V2 (5) 

where "07,^(0 i s the hyper-radial wave function solution of Eq. (|) , ^ [^f^ (O) 
the hyper-spherical harmonics defined in the hyper-sphere of unitary radius, 
and and x the flavor and spin wave function of mixed SU(2) symmetry. 
Let us note that, in order to preserve relativistic covariance, the spin wave 
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functions 



xms = ^= [2 TU - (Tl + IT) T] 



xma = ^= (U - IT) T , (6) 



have to be formulated by means of the appropriate Melosh transformation of 
the ith quark spin wave function, i.e., each individual spin vector must be 
Melosh-rotated: 



Xi = D ' [R M (ki)]Xi = — — r2 ,.,„ xt 



(7) 



where xt are the usual Pauli spinors for the particle i. 

The calculation of the helicity distribution at the hadronic scale /ig can be 
written, according to the expression (uh, as: 



TT M N 



where 



9 (1 -x) 2 

dkj_ Tt(k z ; k | 



B a (m + M NJ ^ 



C a k\ 



-D(k z ,x) (8) 



D(k z ,x) 
k z (x, kj_) 



M Nj ^-k z 



M 



2k 2 



Mn_ 
2 



(kl + m 2 ) l-x 



Ml 



(l-x) 



(9) 
(10) 



and n(k 2 ,Kj_) is the total (unpolarized and flavorless) momentum density 
distribution of the valence quarks in the nucleon calculated making use of 
the eigenfunction -00 o of Eq. (Q) and properly normalized to the number of 
particles. The coefficients B a and C a take the values B u = ~C U — 4 and 
Bd = —Ca = — 1. 

The mass equation (|[) has been solved numerically by fitting the param- 
eters of the potential to reproduce the basic features of the (non strange) 
baryonic spectrum up to ~ 1500 MeV, namely the position o£the Roper res- 
onance and the average value of the l - states. We obtain □: r = 3.3 and 
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Figure 1: Left panel (Fig. la): the polarized distribution Auy(x, /^q) as a function of x: 
the non-relativistic approximation (dot-dashed curve) , and the relativized solution of Eq. (0) 
which neglects Melosh rotations (dotted curve) are compared with the results of a complete 
light-front calculation (full curve). On the right panel (Fig. lb) the distribution Ady (x, /i^) 
(same notation). 

ki = 1.8 fm~ 2 , to be compared with the corresponding non-relativistic i fit 
t = 4.59 and Ki = 1.61 fm~ 2 . The constituent quark masses have been chosen 
m u = m d = m q = M N /3. 

An important outcome of this relativistic treatment is the large amount 
of high momentum components in the wave function, that P|lay an important 
role in the evaluation of transition and elastic form factors u and in the large 
x region in DISH. 

The relevant effects of relativistic covariance are more evident looking at 
the polarized distributions A%(i,/1q) = gi(x,^), AcI v (x,[Iq) = gHx,^) 
where the spin dynamics on the light-front plays a crucial role (Fig. [y). The 
introduction of the Melosh rotations results in a substantial enhancement of the 
responses at large x and in an suppression of the response for 0.1 ~ x ~ 0.5 as 
can be seen from Fig. 1. We show, in the same figure, also the predictions of a 
pure relativized solution obtained by solving numerically Eq. (Q) and neglecting 
the Melosh rotation effects in (||) (dotted curves). Such a calculation retains 
the contribution due to the high momentum components, while the covariance 
requirement on the parton distribution is lost. 

3 Results 

We shall give results in two scenarios characterized by different gluon distri- 
butions (AG) at the hadronic scale: 

i) Scenario A : AG(x, /i 2 ,) = 0. Only quark valence distributions are allowed 
at the hadronic scale. The momentum sum rule determines /i 2 , = 0.094 
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GeV 2 at NLO ([a s (^)/(47r)] NLO = 0.142)0. 



ii) Scenario B: AG(x, //q) = / G(x, /z 2 ,). / is the fraction of polarized gluons 
and has to be considered with the appropriate sign. 

A natural choice for the unpolarized jduon distribution within the present 
approach, has been discussed in refslU'Eil and assumes a valence-like form 
G(x, jLtg) = ^f- \uv{x, /Uq) + d\/(x, /Ltg)] . This definition implies J G(x 7 jj^) dx = 
2 and therefore only 60% of the total momentum is carried by the valence 
quarks at the scale /k 2 ,. In this case /x 2 , = 0.220 GeV 2 at NLO ([a s (^Q)/(47r)]NLO 
= 0.053). If the iduons are fully polarized one has \AG(x, ^)\ — G(x,/j,q). 
Jaffe's suggestionEil in our approximation implies AG(x, /X 2 ,) ~ —0.35 G(x, /x 2 ,) 

In Fig. g the results for the proton structure function g\(x,Q 2 ) = 
\ 12a e 1(9i( x > Q 2 )+9i( x i Q 2 )) are shown and compared with the experimental 
data. The non-relativistic approximation of the present calculation appears to 
reproduce rather poorly the experimental observations. Even the use of non- 
relativistic models which reproduce rather well the baryon spectrum does not 
alter this conclusion, as shown in refO. 

The full covariant calculation leads to theoretical predictions quite close 
to experimental data in the region 0.01 ~ x ~ 0.4 under the assumption of 
a pure valence component at the hadronic scale (scenario A). The calculation 
is parameter- free and the only adjustable parameters (r and ki in Eq. (0)) 
have been fixed to reproduce the low-lying nucleon spectrum as discussed in 
the previous section. The effect of relativistic covariance in the quark wave 
function is mainly associated to the spin dynamics induced by the Melosh 
rotations (|^) □ and these transformations lead to a strong suppression of this 
structure function in the small- a; region (x ~ 0.5). 

In order to introduce the gluons non perturbatively we evolve the unpo- 
larized distributions predicted by the scenario A, up to the scale of scenario 
B where 60% of the total momentum is carried by valence quarks. At that 
scale we substitute the bremsstrahlung gluons by the valence gluons as defined 
previously. Moreover, at that scale the fraction of polarized gluons is chosen 
to be negative according to Jaffe's result Ea. Note that we are maximizing the 
difference with respect to the radiative gluons, because those lead to a positive 
polarization. Fig. 2 leads us to conclude that the low- a: data on g\ do not 
constrain the gluon polarization strongly. If we vary the fraction of polarized 
gluons from 35% to 100% the quality of the agreement is deteriorated in the 

a In fact in ref0 it has been shown that J AG(x, fig) dx < 0. Such inequality does not imply 
AG(x, Hq) < in the whole x-range. 
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Figure 2: The proton polarized structure function at Q 2 = 3 GeV 2 . The full curve represents 
the NLO (MS) results of a complete light-front calculation within a scenario where no gluons 
are present at the hadronic scale (scenario A); the corresponding non-relativistic calculation 
are shown by the dot-dashed line. Scenario B is summarized by the dotted line in the case of 
negative polarized gluon fraction (J AG = —0.7 as discussed in the text), and by thejiashed 

line in the case of positive gluon polarization ( J AG = +0.7). Data are from SMCE3, and 
SLAC(E142)y experiments. 



region 0.01 ~ x ~ 0.4, but only slightly. For larger values of x the valence 
contribution plays a major role and the behavior of the structure functions 
depends largely on the model wave functions. 

Concerning the gluon distribution, our results are summarized in Fig. [|. 
Within scenario A the polarized gluon distributions remain positive at the 
experimental scale as a result of the NLO evolution. On the contrary a large 
amount of negative gluon polarization is predicted within scenario B if one 
assumes anti-aligned gluons at the hadronic scale. The distribution is largely 
dependent on the polarization fraction. The dotted line of Fig. ||a shows 
the consequence, in the deep inelastic regime, of the Jaffe's calculation at low 
energy. The absence of non-perturbative sea polarization results in a huge (and 
probably unrealistic) amount of negative gluon polarization to reproduce the 
neutron data. Such a large amount of gluons is, however, inconsistent with the 
unpolarized gluon distributions as it is shown in Fig. (dotted line again). 

4 Transversity and Relativistic Spin Effects 

In addition to the unpolarized and helicity parton distributions, a third parton 
distribution, the so called transversity h\{x, Q 2 ), is required to get a complete 
twist-2 description of the momentum and spin structure of the nucleonEj. As 
the most remarkable feature hi is chiral-odd and as a consequence it decouples 
from electron-nucleon DIS and other easily accessible hard processes. In doubly 
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Figure 3: Left panel: Polarized gluon distributions at Q 2 = 3 GeV 2 obtained evolv- 
ing at NLO (MS) the polarized partons in both scenarios. Scenario A (full line). Sce- 
nario B is summarized by the dotted line in the case of negative polarized gluon fraction 
(J AG(x, fig) = —0.7 as discussed in the text), and by the dot-dashed line in the case of 
positive gluon polarization ( J AG(x, Hq) = +0.7). Right panel: unpolarized gluon distribu- 
tion at Q 2 = 10 GeV 2 obtained evolving at NLO (DIS) unpolarized partons. Scenario A 
(full line); scenario B (dotted line). For comparison also the results at LO are shown in this 
case (dot-dashed). CTEQ4 NLO (DIS) fit of ref.U: full triangles. 



polarized Drell-Yan processes (pp — > l + l~X) the chirality of the partons that 
annhilate is uncorrelated and, in principle, h\ could be measured (experiments 
are already included in the research program at HERA and RHICE3). 

We will focus on an observable defined as the ratio between double trans- 
verse and double longitudinal asymmetry in the cross section: 



R( n?) = Att 1 = Eg elKjxi, Q 2 M(x 2 ,Q 2 ) + (si ~ x 2 ) 

y !, 2 , v j _ A ^ ^ ^ elgf(x lt Q^(x 2) Q 2 ) + (*i - x 2 ) 

(11) 

where f(9, <j>) is a function of the angles of the final detected lepton, a indicates 
the flavor of the parton, Q 2 is the invariant mass of the produced pair and the 
variables X\ and x 2 are defined in terms of Q 2 , the center of mass energy y/s 
and the quantity y — arcth(Q 3 /Q°): 




(12) 



The physical meaning of h\ is similar to that of g\ but in a transverse 
basis. It is clear that in any non-relativistic description of the nucleon, where 
motion and spin commutes, h\{%,pfy — g1(x,^) at the hadronic scale /io- 
However, Melosh rotation is able to break this degeneracy and h\{x, can be 




Figure 4: Ratio between transverse and longitudinal asymmetries (eq (|ll|)) as a function a) 
of the center of mass rapidity y for Q 2 = 100 GeV 2 and a center of mass energy yfs = 100 
GeV (solid line); b) of the invariant mass of the produced lepton pair (Q 2 ) for a center of 
mass rapidity y = and a center of mass energy s/s = 100 GeV. The dashed line corresponds 
to the case when the Melosh rotation is switched off. The dotted line is the result obtained 
in a fully non-relativistic approach. 



expressed in the same form as gi(x, /j&L eq. (||), but with different coefficients: 
B u = 4, C u = and B d = -1, C d = 00. 

By employing the wave functions calculated in the previous sections and 
NLO evolution we have calculated R(x\, X2, Q 2 ) at a scale Q 2 — 100 GeV 2 
(Fig. |i|a,). Comparing the curves shown in Fig. ||a it is clear that if MR is 
not taken into account, i.e. we start from the equality /i"(x,/Lto) = <7i (^j Mo)) 
the obtained value for R is approximately a factor 2 smaller than when MR 
is considered. The same conclusion holds for a wide kinematic regimen (see 
Fig. ^b). Furthermore, the chosen observable is quite insensitive to the details 
of the mass operator. To check the dependence on the chosen spatial wave 
function we have recalculated R(x\,X2, Q 2 ) with the non-relativistic version 
of the mass equation ^ and results are also shown in Fig. || Therefore we 
can conclude that the chosen observables is a rather direct measure of the 
importance of relativistic spin effects in the nucleon (MR) while it is quite 
insensitive to the details of the spatial wave function. 
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